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We study reaction zones in three different versions of the A + B — > system. For a steady state 
formed by opposing currents of A and B particles we derive scaling behavior via renormalization 
group analysis. By use of a previously developed analogy, these results are extended to the time- 
dependent case of an initially segregated system. We also consider an initially mixed system, which 
forms reaction zones for dimension d < 4. In this case an extension of the steady-state analogy gives 
scaling results characterized by new exponents. 
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The reaction-diffusion process A + B — > has received 
much attention since the work of Toussaint and Wilczek 
0, in which it was demonstrated that density fluctua- 
tions are important when the dimensionality d < 4. At- 
tempts to incorporate the fluctuation effects, including 
the application of renormalization group (RG) methods, 
have proven successful for one-species reactions ^,H,|). 
However, the physics of the two-species reaction is some- 
what different, as the number difference of A and B par- 
ticles is locally conserved. Nonetheless, the main results 
of the RG analysis may be extended to this case |||J . 

In the present work we are concerned with situations 
in which the A and B particles are well segregated, so 
that reactions are then confined to reaction zones, on 
the boundaries between the A- and £>-dominated regions. 
The segregation may be a consequence of the initial con- 
ditions, or arise asymptotically from a mixed initial state 
when d < 4. Both of these cases may be compared to 
that of a steady state formed by directing steady, uni- 
form currents of A and B particles towards each other. 
In all three cases the structures of the reaction zones are 
very similar. 

Reaction zones were first studied by Galfi and Racz |t]] 
in the context of a system with segregated initial con- 
ditions. Their analysis is essentially a mean-field result, 
as we will show, which holds asymptotically for d > 2. 
Later, Cornell and Droz treated the case of d < 2 by 
using scaling arguments motivated by a renormalization 
group picture g]. Also, reaction zones in initially ho- 
mogeneous systems were studied by Leyvraz and Redner 
P,[To[. In this paper we use exact RG methods both to 
confirm previous results and to derive new ones for the 
exponents which characterize the reaction zones. 

The model we shall use for the A + B — » reaction 
is one of particles undergoing continuous-time random 
walks on a hypercubic lattice. We consider general, non- 
zero values of the diffusion constants Da and Db- If an 
A and a B particle are together on a lattice site, then 
they annihilate with some characteristic reaction rate A. 
In the field-theoretic approach it is convenient to allow 



multiple occupancy of lattice sites, regardless of parti- 
cle type, but this is not expected to alter the universal 
properties. 

First we consider the steady state reached when equal 
currents J of A and B particles are directed towards each 
other. The average densities then vary only in the direc- 
tion of the currents. In this case the width of the reaction 
zone goes as a power of J, as J — » 0, 



d<2 



j-l/3 



d>2. 



(1) 



The typical nearest-neighbor distance in the reaction 
zone £ TZ scales as 



4z 



j-i/(d+i) d < 2 

j-2/3d d > ^ 



(2) 



The second example we consider is that of initially seg- 
regated particles. Consider a system prepared with only 
A particles in the region x < and only B particles for 
x > 0. The behavior in the reaction zone again scales, 
now with respect to time, with w ~ t a and £ TZ ~ V ' . The 
densities in the reaction zone, where \x — x c \ < w, have 
the scaling form 



(a),(b) ~t-i d F a 



(3) 



where x c oc t 1 / 2 is the reaction zone center. The angular 
brackets will be defined precisely below. The reaction 
rate R also has the scaling form 



R(x,t) ~ t~ p G 



(4) 



The values of the exponents are given in table |. 

The final example is that of random, homogeneous ini- 
tial conditions, with equal densities no of A and B par- 
ticles. It has been shown in the mathematical treatment 
of Bramson and Lebowitz [fy] that this system exhibits 
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TABLE I. The exponents, as defined in the text, for the 
scaling behavior in the cases of segregated and homogeneous 
initial conditions. 
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asymptotic segregation for d < 4. We present a sim- 
ple derivation of this result from an effective field theory, 
which is valid for 2 < d < 4. We argue that reaction 
zones which result from the segregation can be studied 
in the same fashion as in the previous problems. As a re- 
sult, we predict the same scaling forms as given by Eqs. 
(|) and (§), with new values for the exponents, which are 
listed in rows 3 and 4 of table |. 

The characteristic nearest-neighbor distance in the re- 
action zone £ vz is found from the scaling form of the den- 
sities in the reaction zone ([|) . Consider the total number 
of particles N in the volume of the reaction zone for a 
given interfacial surface area Sd—i, 



dx (a(x) + b(x)) 



S d ^t a -< d . (5) 



Since the volume of this region is Sd-it a , then the aver- 
age volume per particle goes as V d . If one assumes that 
there is just one length scale describing typical nearest- 
neighbor distance, both in the direction of the interface 
and perpendicular to it, and also for both like and unlike 
particles, then this length is given by £ vz ~ t 1 . A conve- 
nient definition of £ rz , which follows from the derivation 
above, is £ TZ = {a{x c ))~ 1 / d . 

The model for A + B — > can be mapped to a field 
theory, which is useful for the application of RG meth- 
ods §|TJ. We summarize here the results of a study of 
this field theory which pertain to the problem at hand. 
A more detailed account will appear in |(|. One fea- 
ture, which is general to all irreversible reaction-diffusion 
systems, is that there are no diagrams which dress the 
propagator. This means that there is no wavefunction 
renormalization, and consequently no anomalous dimen- 
sion for the fields. The significance of this will be noted 
below. 

There is only one coupling constant Ao in the field the- 
ory (in the notation of [g| ) , which is given in terms of the 
original microscopic parameters as Ao = Xh d /D, where 
D = (Da + Db)/2 and h is the lattice size. The coupling 
Ao is found to be irrelevant, in the RG sense, for d > 2, 
marginal for d — 2, and relevant for d < 2. 

In general one can derive from the field theory the 
equations of motion 



dt{a) =D A V 2 (a) -R d t (b) = D B V 2 (b) - R (6) 

where R — \$D{ab) is the reaction rate. The angular 
brackets refer to averages over the stochastic processes of 
diffusion and reaction, but not over the initial conditions. 
Averages over the initial conditions will be denoted with 
a bar: (a). Note that the quantity (a — b) obeys a simple 
diffusion equation when Da = Db- 

When the d > 2 the irrelevance of the coupling leads 
to the asymptotic result that R ~ T(a)(b). Then the 
equations of motion (||) are simply differential equations 
with rate constant T, which is the starting point in the 
analysis of Galfi and Racz . Therefore their results are 
applicable for d > 2. When d < 2 one finds that the 
coupling Ao can be renormalizcd exactly, allowing us to 
derive the scaling behavior of the reaction zones. 

Consider the steady state reached by imposing currents 
J a = Jx at x — —L and = — Jx at x = L, where 
L 3> w. From the analysis in M one finds for d > 2 that 



R ~ (a)(b) ^ J 4 / 3 /M 1/3 ), 



(7) 



from which the exponents for d > 2 given in Eqs. (Q) and 
(||) follow directly. 

For d < 2 one must consider the effects of renor- 
malization. We work in units of time where the aver- 
age diffusion constant D = 1, and define the parameter 
6 = (Da — Db)/(Da + Db)- One can show that 8 is not 
renormalized, for the same reason that there is no wave- 
function renormalization. In these units the dimension 
of the current is [J] = (length) -1- ' 1 . The width w can 
depend in the steady state on Ao, J, and S only. 

In the field-theoretic RG, the Ao-dependence is traded 
for a dependence on a renormalized coupling A/j, defined 
at an arbitrary length scale J 1 ^ d+1 \ We then consider 
w as depending on J, 6, Jq, and the the dimensionless 
renormalized coupling gn 
satisfies an RG equation whose solution is 



XRj ( d -2)/( d +l)_ ^ 



w(J,g r ,J ,5) 



-i/(d+i) 



<(jo,g(J/Jo,ga),6) (8) 



where, as J — > 0, the running coupling g — > g* = 0(2— d), 
independently of the initial value of A . When d = 2 the 
running coupling goes as g ~ B/hi\J\ for small J. The 
simple form of the overall scale factor is a result of the 
absence of wave-function renormalization. 

The significance of Eq. (||) is that, in the limit of small 
J, the only J dependence on the right-hand side comes 
from the overall scale factor. Therefore any length will 
scale asymptotically as w ~ £ rz ~ J -1 /( rf + 1 ) j giving the 
results stated Eqs. @) and (§) for d < 2. When d = 2 the 
running coupling g is J dependent. However, it can be 
shown that the right-hand side is finite in the limit g — -> 
Q. Therefore the leading term for small J is unchanged, 
and there are no logarithmic corrections. It should be 
noted that while this result holds for general \S\ < 1, the 
limit of | S | — » 1 does not necessarily commute with the 
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FIG. 1. The profile of a reaction zone. 



asymptotic limit in our analysis, so our results are only 
for the case when both species are mobile. 

The asymptotic limit of any dimensionful quantity can 
be found by the same technique, using only dimensional 
analysis. Therefore the densities in the reaction zone 
have the scaling form 



(a(x)),(b(x)) ~ J^^Fa^xJ 1 /^) 



and the reaction rate 



R(x) ~ J^+W+^GixJ 1 '^). 



(9) 



(10) 



The scaling functions F a ,t and G are expected to be 
universal, up to metric factors, and should be calcula- 
ble within a (2 — d)-expansion. Note that R oc (ab) ~ 
j(2-d)/(d+i)^^ for gmall j Thig regult can be demon _ 

strated explicitly within the original field theory Q . 

An equivalent and less formal way to state these re- 
sults is that there are only two input parameters with 
dimension: J, and Ao- However, for d < 2, in the asymp- 
totic limit the theory is independent of the initial value 
of Ao, so the remaining J dependence can be determined 
on dimensional grounds. 

These results were found by Cornell and Droz ||, by 
use of an RG- motivated scaling analysis. Our work pro- 
vides a quantitative confirmation of the lack of anomalous 
dimension and asymptotic lack of Ao dependence. They 
performed simulations for d = 1,2,3 which are in good 
agreement with the predictions. 

We now turn to the problem of segregated initial con- 
ditions. At t = the boundary is at x — 0, and the 
A and B particles are randomly distributed within their 
regions, with initial densities ao and bo respectively. A 
profile of the densities at a later time t is sketched in Fig. 
|l|. The densities are depleted out to a range Wd ~ t 1 / 2 , 
which is the length over which particles will have had 
a chance to diffuse into reaction range. Provided that 
a < 1/2, which will be verified self-consistently, then 
asymptotically w <C Wd- Consequently one finds in the 
depletion regions w -C \x — x c \ -C Wd that R « and 
Eq. (^) reduces to the diffusion equation. As a result, it 
follows that the density profiles in this depletion region 
are linear in x — x c . The slope can be determined by ob- 
serving that the A particle density goes from ao to zero 
in a range Wd, so the slope — ao/Wd ~ —t~ x / 2 . 

Implicit in the analysis of Galfi and Racz |j is the 
observation that this linear depletion regime provides the 
same boundary conditions for the reaction zone as was 



found in the steady-state system. In that case the slope 
in the linear region (w <C \x\ < L) is given by d x {a) = 
—J. Therefore we can apply the results obtained from 
the steady state to the time-dependent reaction zone by 
making the scaling substitution J ~ t^ 1 ! 2 . The result 
is the scaling forms given by Eqs. (|) and Eq. (|) with 
the exponents shown in table Q. The exponent 7 for the 
nearest-neighbor distance £ rz is, to our knowledge, a new 
result. 

The exponents a and /3, which describe w and R, were 
derived in the analysis of Cornell and Droz Q . They find 
reasonable agreement with these predictions in d = 1 
simulations, although the asymptotic region is difficult 
to obtain. Other simulations in d = 1 find evidence for 
multi-scaling |]l3| . Since the parameters in these simula- 
tions are at the opposite extreme from the weak coupling 
expansion implicit in the field-theoretic approach, it is 
conceivable that they might fall into a separate univer- 
sality class. However, more extensive simulations seem to 
indicate that ordinary scaling is recovered asymptotically 

The final case we consider is the system with homoge- 
neous initial conditions. The A and B particles are ran- 
domly distributed with equal initial densities no. Start- 
ing from Eq. (Q) one can show that the average over ini- 
tial conditions yields a density (a) ~ t~ d / i for d < 4 
H . As argued by Bramson and Lebowitz, the A and B 
particles segregate asymptotically, which may be shown 
directly from Eq. (||) for 2 < d < 4 as follows. We use 
the notation (a) — > a, since R = Tab for d > 2. 

Given that / 2 < / 2 , where / is a real quantity, then 

1- 



min(a,6) 2 < - (a 2 + b 2 - (a + b)\a - b\) . (11) 

Since a and b correspond to the physical density for 
a particular initial condition, they are everywhere non- 
negative. Therefore a + b > \a — b\ at every point (x, t), 
which implies that 

1 



min(a, 6) 2 < - (a 2 + b 2 - (a - b) 2 ) = ab. (12) 



From Eq. 



it follows that 
1 



ab - 



d t a = 0(t- 1 - d ^ i ), 



so one has the result 



min(a,6) < O^- 1 / 2 ^/ 8 ) < t~ d / 4 



(13) 



(14) 



when d < 4. Therefore the local minority density goes to 
zero relative to that of the majority species. A sketch of 
the segregated domains in d = 2 is shown in Fig. g. 

In this case we may make a similar analogy to the 
steady-state case, except now the height at the edge of 
the depletion region is time-dependent. The rectangle in 



Fig. 
will 



gives the cross-section of a reaction zone, which 
lave same form as Fig. |[ The time dependence of 
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FIG. 2. A sketch of the percolated domain configuration in 
d = 2. The domains have a characteristic length scale which 
grows in time as t 1 ^ 2 . 



the height is in general complicated, but we argue that 
in most regions it will scale with the bulk density. The 
picture one has is that the densities in the bulk will be 
fairly uniform since there the particles diffuse without 



-d/4 



SO 



reacting. This uniform density must scale as t 
that averages over the reaction zones of the whole sys- 
tem will yield exponents which can be derived simply by 
taking ao,&o — * t~ d l 4 . The resulting steady state analog 
is given by J ~ t-( d + 2 )/ 4 . This gives the scaling results 
Eqs. (||) and (||), with the exponents listed in table | for 
the homogeneous case. 

The exponents a and f3 are, to our knowledge, new re- 
sults. We note that a — * 1/2 from below as d — > 4. This 
is consistent with the view that segregation is breaking 
down, since the reaction zones are then on the same scale 
as the domains. The exponent 7 was derived by Leyvraz 
and Redner |],|l(j by a different method. They were in- 
terested in the quantity £ab, which is the characteristic 
nearest-neighbor distance when the nearest neighbor of 
a particle is of the opposite species. This quantity can 
be measured directly in simulations. Since the particles 
are segregated, the contributing particles to the distri- 
bution for £ab must be within the reaction zone, and 
therefore £ab ~ ^rz ~ t~> '. Leyvraz and Redner found 
from both scaling arguments and simulations values of 
7 which agree with our predictions for d — 1,2. How- 
ever, they argue for d > 2 that 7 = 1/4, and show the 
results of d — 3 simulations. Our prediction for d = 3 is 
7 = 5/18, which is very close to t 1 / 4 . It appears from the 
data of that either value of 7 is an equally good fit. 
A numerical value for a fit to the data is not provided. 
The significance of this discrepancy is that the nearest- 
neighbor distance in the bulk (a) ^ ~ i 1 / 4 for d < 4. 
We claim that £ IZ represents a new length scale for all 
d < 4. 

A check on the above results is to calculate the charac- 
teristic domain size. Integrating the equations of motion 
(p|) over the entire system of size V gives 



d d xd t {a) 



Vr x - d ^ ~aJ dx x R{x±) ~ At 



-0+a 



(16) 



where A is the interfacial area of the domain boundaries. 
This first result follows from the fact that the only sig- 
nificant contribution to the reaction rate comes from the 
reaction zones, the second from Eq. (Q). This leads to 
the characteristic domain size 



VIA ~ t~ (3+a+1+d/i - t 1/2 



(17) 



for all d < 4. This result, while often assumed, follows 
for d < 2 only because R ^ (a) (b) . 
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d d x(ab) 
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or equivalently 
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